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ABSTRACT 

A  post-CME  current  sheet  (CS)  is  a  common  feature  developed  behind  an 
erupting  flux  rope  in  CME  models.  Observationally,  white  light  observations 
have  recorded  many  occurrences  of  a  thin  ray  appearing  behind  a  CME  eruption 
that  closely  resembles  a  post-CME  CS  in  its  spatial  correspondence  and  mor¬ 
phology.  UV  and  X-ray  observations  further  strengthen  this  interpretation  by 
the  observations  of  high  temperature  emission  at  locations  consistent  with  model 
predictions.  The  next  question  then  becomes  whether  the  properties  inside  a 
post-CME  current  sheet  predicted  by  a  model  agree  with  observed  properties.  In 
this  work,  we  assume  that  the  post-CME  CS  is  a  consequence  of  Petschek-like 
reconnection  and  that  the  observed  ray-like  structure  is  bounded  by  a  pair  of  slow 
mode  shocks  developed  from  the  reconnection  site.  We  perform  time-dependent 
ionization  calculations  and  model  the  UV  line  emission.  We  End  that  such  a 
model  is  consistent  with  SOHO/UVCS  observations  of  the  post-CME  CS.  The 
change  of  Fe  XVIII  emission  in  one  event  implies  an  inflow  speed  of  ~10  km/s 
and  a  corresponding  reconnection  rate  of  Ma  ~  0.01.  We  calculate  the  expected 
X-ray  emission  for  comparison  with  X-ray  observations  by  Hinode/XRT,  as  well 
as  the  ionic  charge  states  as  would  be  measured  in-situ  at  1  AU.  We  find  that  the 
predicted  count  rate  for  Hinode/XRT  agree  with  what  was  observed  in  a  post- 
CME  CS  on  April  9,  2008,  and  the  predicted  ionic  charge  states  are  consistent 
with  high  ionization  states  commonly  measured  in  the  interplanetary  CMEs.  The 
model  results  depend  strongly  on  the  physical  parameters  in  the  ambient  corona, 
namely  the  coronal  magnetic  held,  the  electron  density  and  temperature  during 
the  CME  event.  It  is  crucial  to  obtain  these  ambient  coronal  parameters  and  as 
many  facets  of  the  CS  properties  as  possible  by  observational  means  so  that  the 
post-CME  current  sheet  models  can  be  scrutinized  more  effectively. 
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1.  Introduction 

Most  models  for  coronal  mass  ejections  (CMEs),  regardless  of  what  initiates  the  CME, 
predict  a  current  sheet  (CS)  that  develops  beneath  an  erupting  flux  rope  due  to  the  stretching 
of  the  overlying  coronal  held  (e.g.,  Lin  &  Forbes  2000,  Linker  et  al.  2003,  Lynch  et  al.  2004, 
MacNeice  et  ah  2004,  Manchester  et  al.  2004).  Magnetic  reconnection  in  the  CS  reduces 
tension  restraining  the  outgoing  flux  rope  and  at  the  same  time  produces  post-CME  loops 
beneath  the  reconnection  point.  Outflows  along  the  CS  and  a  temperature  higher  than 
that  in  the  ambient  corona  are  expected  inside  the  CS  as  magnetic  energy  is  converted  to 
kinetic  and  thermal  energy  due  to  reconnection.  This  standard  hare-CME  picture  (Figure 
la)  is  supported  by  observations  such  as  loop-top  hard  X-ray  sources  (e.g.  Masuda  et  al. 
1994),  upward  growth  of  hare  loops  (e.g.  Svestka  1996),  separating  hare  ribbons,  and  hotter 
post-hare  loops  lying  higher  than  the  cooler  post-hare  loops  (e.g.  van  Driel-Gesztelyi  et 
al.  1997).  However,  few  observational  signatures  of  the  CS  above  or  near  the  X-point 
were  reported  until  recently.  Sui  et  al.  (2003,2004)  reported  signatures  of  a  CS  at  both 
sides  of  the  reconnection  site  from  X-ray  observations  by  the  Reuven  Ramaty  High  Energy 
Solar  Spectroscopy  Imager  (RHESSI).  Savage  et  al.  (2010)  reported  a  post-CME  CS  feature 
observed  by  Hinode/XRT  with  downhows/uphows  that  allow  the  X-point  to  be  located. 
Innes  et  al.  (2003b)  and  Wang  et  al.  (2007)  reported  high  speed  outflows  in  opposite 
directions  away  from  the  reconnection  site  based  on  Doppler  shift  signatures  of  Fe  XIX  and 
Fe  XXI  lines  observed  by  Solar  Ultraviolet  Measurements  of  Emitted  Radiation  (SUMER) 
on  the  Solar  and  Heliospheric  Observatory  (SOHO).  White  light  (WL)  observations  of  CME 
events  by  Solar  Maximum  Mission  (SMM)  and  Large  Angle  Spectrometric  Coronagraph 
(LASCO)  on  SOHO  have  recorded  many  occurrences  of  a  thin  ray  appearing  behind  a  CME 
eruption  (we  will  call  it  ‘WL  ray’  hereafter)  that  closely  resembles  a  post-CME  CS  in  its 
spatial  correspondence  and  morphology  (e.g.  Webb  et  al.  2003,  Ko  et  al.  2003,  Lin  et 
al.  2005).  Observations  by  Ultraviolet  Coronagraph  Spectrometer  (UVCS)  on  SOHO  of 
high-temperature  emission  (3-6  million  degrees)  that  lies  along  such  WL  ray  in  between  the 
erupting  CME  and  post-CME  loops  provide  strong  support  for  the  CS  interpretation,  as 
opposed  to  being  just  a  usual  streamer  seen  edge-on  (Ciaravella  et  al.  2002,  Ko  et  al.  2003, 
Bemporad  et  al.  2006,  Lee  et  al.  2006,  Ciaravella  &  Raymond  2008).  Figure  2  shows  two 
examples  of  such  events. 

Even  though  there  is  strong  observational  support  for  the  observed  WL  ray  being  the 
post-CME  CS,  that  interpretation  remains  under  scrutiny.  One  major  cause  of  skepticism 
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arises  from  the  observed  thickness  of  the  WL  ray,  as  well  as  the  spatial  extent  of  the  high- 
temperature  emission  observed  by  UVCS,  which  are  much  larger  than  what  most  recon¬ 
nection  models  predict.  The  observed  thickness  is  of  the  order  of  104  —  105  km  (Lin  et  ah 
2007,2009,  Vrsnak  et  al.  2009)  which  is  orders  of  magnitude  larger  than  theoretical  esti¬ 
mates  of  the  thickness  for  the  diffusion  region  and  Sweet-Parker  CS  (see  Bemporad  2008, 
Lin  et  al.  2009)  even  for  the  largest  estimates  of  anomalous  resistivity  currently  available, 
and  taking  projection  effects  into  account.  Bemporad  (2008)  and  Lin  et  al.  (2009)  dis¬ 
cussed  several  reconnection  schemes  with  turbulence  or  instability  (see  references  therein) 
that  may  produce  CS  thickness  comparable  to  what  was  observed,  such  as  turbulence 
and  stochastic/fractal/time-dependent  Petcheck-type  reconnection.  Model  calculations  have 
shown  that  the  effect  of  thermal  conduction  can  create  a  thermal  ‘halo’  around  the  current 
layer  (Yokoyama  &  Shibata  1997,  2001;  Seaton  &  Forbes  2009).  This  would  contribute  to 
the  extent  of  the  high-temperature  emission  observed  by  UVCS,  although  not  for  the  WL 
emission  (which  depends  mainly  on  the  electron  density).  The  electric  potential  of  the  slow 
mode  shocks  might  inhibit  transport  of  energy  into  the  upstream  flow,  reducing  this  effect. 
To  see  which  theories/models  are  more  consistent  with  or  to  reject  the  CS  interpretation, 
the  theories  or  models  for  the  CS  must  predict  some  physical  quantities  that  can  be  tested 
by  these  observations. 

We  should  not  forget  that  there  is  another  piece  of  observational  evidence  that  can 
be  used  to  scrutinize  the  CS  interpretation  and  constrain  viable  theories/models.  As  men¬ 
tioned  above,  SOHO/UVCS  has  observed  several  CME  events  exhibiting  high-temperature 
emission  from  the  [Fe  XVIII]  A974  line  (formation  temperature  of  6  million  degrees)  at  he¬ 
liocentric  heights  of  1.5-1. 7  R0  that  lay  along  the  line  connecting  the  erupting  CME  and  the 
associated  post-CME  loops  (e.g.  see  examples  in  Fig. 2).  Its  location  and  the  timing  of  its 
brightening  following  the  eruption  strongly  favor  the  interpretation  of  a  post-CME  CS.  It 
is  then  important  to  go  one  step  beyond  morphology  and  find  out  if  a  given  CS  model  can 
predict  physical  quantities,  such  as  density  and  line  intensities,  that  agree  with  the  obser¬ 
vations.  And  if  so,  how  the  observed  quantities  would  imply  the  physical  conditions  for  the 
associated  reconnection  within  the  context  of  a  given  model. 

It  is  important  to  note  that  the  analysis  of  emission  lines  from  current  sheets  has  so  far 
assumed  that  the  plasma  is  in  ionization  equilibrium,  but  this  can  be  a  poor  assumption. 
For  instance,  the  ionization  time  scale  for  Fe  XVIII  is  on  the  order  of  1010/ne  seconds,  or 
nearly  1000  seconds  for  the  densities  estimated  from  UVCS  observations,  while  the  flow  time 
in  the  current  sheet  should  be  on  the  order  of  a  few  tenths  of  a  solar  radius  divided  by  an 
ouflow  speed  on  the  order  of  1000  km  s-1,  which  can  be  several  times  smaller.  It  is  therefore 
important  to  consider  time-dependent  ionization  when  interpreting  emission  line  intensities. 
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Vrsnak  et  al.  (2009)  proposed  a  working  hypothesis  for  the  observed  post-CME  CS 
(i.e.  WL  ray),  and  provided  the  Erst  attempt  to  calculate  the  electron  density  within  the 
CS,  and  compare  with  the  observations.  This  model  assumes  a  steady  state  Petschek-like 
reconnection  scheme  (Petschek  1964).  The  reconnection  takes  place  in  the  diffusion  region 
(DR)  and  a  pair  of  standing  slow-mode  shocks  (SMSs)  develop  and  extend  out  from  the 
DR.  At  the  SMS  crossing,  the  inflowing  coronal  material  is  compressed  and  heated  in  a  way 
that  is  mainly  determined  by  the  external  plasma  f3,  and  the  velocity  of  the  outflowing  jet 
within  the  region  bounded  by  the  SMSs  is  approximately  equal  to  the  external  Alfven  speed 
(Aurass  et  al.  2002).  They  found  that  the  observed  electron  density  and  morphology  of  the 
post-CME  CS  (including  thickness)  can  be  successfully  explained  by  their  model.  In  this 
paper,  we  adopt  their  model  and  calculate,  for  the  first  time,  the  expected  UV  and  X-ray 
emission  inside  the  CS.  Strictly  speaking,  within  the  context  of  this  model,  the  so-called 
post-CME  CS  seen  as  a  ray-like  structure  in  white  light  observations  is  actually  the  region 
bounded  by  the  SMSs.  This  region  is  not  a  ’current  sheet’  per  se,  and  the  electric  currents 
are  concentrated  only  in  the  DR  and  SMSs  (Vrsnak  et  al.  2009).  Note  that,  however,  we 
will  use  ‘post-CME  CS’,  ‘CS’,  or  ‘WL  ray’  interchangeably  throughout  this  paper  to  stand 
for  this  region. 

Section  2  describes  the  models  in  detail.  We  adopt  two  approaches.  One  is  the  ‘Fully- 
mixed’  model  which  is  the  direct  adaptation  of  the  model  laid  out  in  Vrsnak  et  al.  (2009). 
The  other  is  the  ‘Streamline  model’  which  assumes  that  the  outflow  jets  out  of  the  SMS 
crossings  at  different  locations  do  not  mix  with  each  other.  We  describe  the  methods  for 
calculating  the  electron  density,  electron  temperature  and  time-dependent  ionic  fractions 
within  the  CS  for  both  models.  Section  3  presents  the  predicted  UV  and  X-ray  emission, 
and  the  expected  frozen-in  charge  states  as  would  be  measured  in-situ  at  1  ALL  We  discuss 
the  results  in  Section  4. 


2.  Modeling  the  Post-CME  Current  Sheet 

Figure  lb  illustrates  the  idea  of  the  model  of  Vrsnak  et  al.  (2009)  that  is  based  on 
the  steady-state  Petschek  reconnection  scheme.  As  the  eruption  stretches  the  oppositely 
directed  coronal  held  beneath  the  flux  rope,  reconnection  takes  place  at  the  DR  under  suitable 
conditions.  A  pair  of  standing  SMSs  then  develops  out  of  the  DR  both  above  and  below.  In 
this  paper,  we  will  concentrate  on  the  region  bounded  by  the  pair  of  SMSs  above  the  DR  in 
this  vertical  CS  configuration.  Coronal  material  that  crosses  the  SMS  is  compressed,  heated, 
accelerated,  and  forms  an  outflow  jet  above  the  DR.  The  extent  of  compression  and  heating 
depends  on  the  ambient  (i.e.  coronal)  plasma  f3  at  the  location  of  the  SMS  crossing  as  in 
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Aurass  et  al.  (2002): 

n2  _  5(1  +  (3) 
rn  ~  2  +  5  f3 

and 


(1) 


T2  2 

—  =  H - 

T\  ^5/3 


(2) 


where  ni,Ti  and  n2,T2  are  the  electron  density  and  electron  temperature  in  the  ambient 
corona  and  after  the  SMS  crossing,  respectively.  (3  is  the  ratio  of  gas  to  magnetic  pressure 
in  the  ambient  corona  for  fully  ionized  plasma  (Mann  et  ah  1999): 


l.92riikTi 

(Bcor /870 


(3) 


where  Bcor  is  the  ambient  coronal  magnetic  held  at  the  SMS  crossing.  The  outflow  speed 
out  of  the  SMS  crossing  is  approximated  to  be  the  Alfven  speed  for  the  inflowing  material: 
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To  calculate  the  physical  properties  inside  this  post-CME  CS,  one  needs  to  take  into 
account  all  outflows  from  the  SMS  crossings  along  the  SMS  (Fig. lb).  At  this  point,  we 
will  investigate  two  models.  One  is  the  ’Fully-mixed’  model  (‘Model  1’,  Fig.lc)  in  which 
plasma  that  just  comes  inside  the  CS  through  SMS  mixes  thoroughly  with  material  that 
comes  from  below,  as  in  Vrsnak  et  al.  (2009).  This  is  the  ‘fluid’  approach  and  the  electron 
conductivity  is  large  across  the  CS  so  that  thermal  equilibrium  is  achieved  between  the 
incoming  and  outflow  material  at  a  given  height.  The  other  is  an  extreme  opposite  of  the 
Erst  model,  the  ’streamline’  model  (‘Model  2’,  cp.  Fig. lb),  in  which  plasma  that  just  comes 
inside  the  CS  through  SMS  forms  its  own  ‘streamline’  and  does  not  mix  with  the  material 
that  comes  in  at  other  SMS  crossings.  This  would  be  the  case  if  certain  conditions  inside 
the  CS,  e.g.  turbulence,  prevent  these  flows  (i.e.  streamlines)  from  efficiently  changing  the 
thermodynamical  properties  of  each  other.  We  discuss  the  two  models  in  more  detail  below. 

Both  models  adopt  the  same  ambient  coronal  conditions.  As  shown  above,  the  properties 
inside  the  CS  depend  on  the  ambient  plasma  (3,  which  in  turn  depends  on  the  electron 
density,  electron  temperature  and  magnetic  held  profiles  in  the  ambient  corona.  For  the 
coronal  magnetic  held,  we  adopt  the  empirical  model  by  Dulk  &  McLean  (1978): 

Bcor(r)  =  0.5(r  —  1)-L5G  (5) 

where  r  is  the  heliocentric  distance  from  the  Sun  in  R0.  Strictly  speaking,  this  formula  is 
applicable  for  approximately  1.02  <  r  <  10  R0  but  we  will  use  this  formula  up  to  20  R0  (see 


6 


also  Vrsnak  et  al.  2002),  the  upper  boundary  of  our  calculations.  The  accuracy  of  this 
formula  between  10  and  20  R0  is  not  a  concern  for  our  study  since,  as  we  will  see  below,  1) 
the  emission  calculated  here  only  have  observational  data  available  at  locations  much  lower 
than  10  R0  (because  the  emission  drops  with  n2),  and  2)  the  change  of  time-dependent  ionic 
fractions  occur  mostly  below  10  R0.  This  formula  gives  a  value  of  1.6  x  10  4  Gauss  (16 
nT)  at  1  AU  which  is  also  reasonable  (e.g.  see  ACE  L3  summary  plots  at  the  ACE  Science 
Center,  http://www.srl.caltech.edu/ACE/ASC/index.html).  To  investigate  the  dependence 
on  different  coronal  electron  density  and  temperature  profiles,  we  use  two  ncor  and  two 
Tcor  profiles.  For  the  coronal  electron  density  profiles,  one  (denoted  as  ‘Nl’)  is  a  ’hybrid’ 
profile  which,  for  r  <  5.66  R0,  is  obtained  for  a  streamer  during  the  SPARTAN  201-1 
mission  (Guhathakurta  &  Fisher  1995).  Beyond  5.66  R0,  we  linearly  interpolate  between 
the  measured  values  at  5.66  R0  and  a  value  of  211.2  cm-3  at  43  R0  (0.2  AU)  which  is  based 
on  Mann  et  al.  (1999): 

nCOr{r )  =  6.53  J  cm"3.  (6) 

for  0.2  AU  <  r  <  5  AU.  The  second  ne(r )  (denoted  as  ‘N2’)  is  adopted  from  LeBlanc  et  al. 
(1998)  which  is  of  an  analytic  form  of 

ncor(r)  =  3.3  x  105r“2  +  4.1  x  106r“4  +  8.0  x  107r"6  cm"3.  (7) 


Note  that  both  profiles  are  based  on  observations.  For  the  coronal  electron  temperature 
profiles,  we  assume  a  general  form  (in  K): 


J  T0  for  1  R0<  r  <  2  R0; 

\T0/(0.5r)a92  for  r  >  2  R0. 


(8) 


The  index  0.92  is  based  on  Totten  &  Freeman  (1995)  for  a  spherically  expanding  wind.  For 
T0,  we  will  compare  results  for  T0  =  106  K  (‘Tl’)  and  2  x  106  K  (‘T2’).  The  solar  wind  speed 
profile  is  then  taken  as: 

,  x  „  1  / 215\  2 

vsw{r )  =  2  x  108 - ( -  cm/s.  (9) 

nCor{r)  V  r  J 

based  on  the  continuity  equation  nvr2  =constant  with  nv  =  2  x  10s  cm"2s_1  at  1  AU. 

Figure  3  plots  the  adopted  Bcor,  ncor,  and  Tcor  profiles,  and  the  corresponding  plasma 
/ 3 .  Also  plotted  (‘lMK_Parker’)  are  those  adopted  by  Vrsnak  et  al.  (2009)  using  a  Parker 
wind  model  (Mann  et  al.  1999)  with  1  MK  isothermal  corona,  and  the  same  Bcor{r )  (Dulk 
&  McLean  1978).  Note  that  none  of  the  actual  coronal  parameters  in  the  immediate  vicin¬ 
ity  of  the  SMSs  are  known  for  certain,  and  the  coronal  conditions  are  likely  to  vary  from 
one  event  to  another.  Furthermore,  Vrsnak  et  al.  (2009)  found  that  the  coronal  density  is 
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depleted  from  the  pre-CME  corona  in  the  vicinity  of  the  WL  ray.  In  any  case,  we  believe 
that  these  coronal  profiles  adopted  here  are  reasonable  and  can  be  taken  as  generic  choices. 
Nonetheless,  it  is  necessary  to  investigate  the  effect  of  different  coronal  parameters  on  the 
calculated  CS  properties.  Therefore  we  calculate  the  model  results  for  three  cases:  1)  Nl+Tl 
(denoted  as  ‘N1T1’),  2)  N1+T2  (denoted  as  ‘N1T2’),  and  3)  N2+T2  (denoted  as  ‘N2T2’) 
while  keeping  the  same  Bcor(r).  Cases  ‘N1T1’  and  ‘N1T2’  are  used  to  investigate  effects 
from  different  coronal  temperature  profiles.  Cases  ‘N1T2’  and  ‘N2T2’  are  used  to  investigate 
effects  from  different  coronal  density  profiles.  We  take  N1T2  as  the  reference  case  for  com¬ 
parison  among  different  cases  and  models.  The  lower  right  panel  of  Fig. 3  plots  the  plasma 
/?,  which  determines  the  jump  condition  across  the  SMS  (Eqs.(l)  and  (2)),  for  these  3  cases. 
The  (5  in  Vrsnak  et  al.  (2009,  ‘lMK_Parker’)  is  also  plotted  for  comparison. 


2.1.  Model  1:  The  Fully-mixed  Model 

Model  1  (Fig.lc)  follows  the  same  thermodynamic  calculations  as  Vrsnak  et  al.  (2009), 
except  that  we  calculate  the  CS  properties  with  three  cases  of  coronal  parameters  described 
above.  Readers  are  referred  to  the  Appendix  in  Vrsnak  et  al.  (2009)  for  a  complete  descrip¬ 
tion  of  the  model.  We  only  briefly  describe  the  concept  here.  For  a  given  cell  (A R  in  Fig.lc) 
at  a  given  height  Rt  above  the  DR,  plasma  comes  into  the  CS  through  the  SMS  with  density 
n2  (Eq.(l))  into  a  given  volume  (Vsms)  that  is  determined  by  the  ambient  coronal  param¬ 
eters.  There  is  also  material  coming  from  the  cell  below  into  a  volume  ( Vthru )  determined 
by  spherical  expansion.  In  this  model,  the  outflow  (i.e.  exhaust  flow)  speed  is  taken  as  the 
local  value  of  va  V  ^sw  (for  details  see  Vrsnak  et  al.  2009)  but  vsw  is  almost  always  negligible 
in  the  considered  height  range.  Under  the  assumption  that  the  ‘incoming’  and  ‘through’ 
material  thoroughly  mix  with  each  other,  the  electron  density  flowing  out  of  this  cell  into 
the  cell  above  it  can  be  calculated  from  mass  conservation,  along  with  the  CS  geometry.  The 
electron  temperature  in  the  cell  is  calculated  to  be  the  average  of  the  temperature  after  the 
SMS  crossing  (i.e.  T2,  Eq.(2))  and  that  coming  from  the  cell  below  under  adiabatic  cool¬ 
ing,  weighted  by  the  mass  inside  Vsms  and  Vthru  respectively.  The  density  and  temperature 
profiles  along  the  CS  (above  the  DR)  can  then  be  calculated  iteratively,  and  they  depend 
on  the  ambient  coronal  parameters  and  the  location  of  the  DR.  Since  the  magnetic  field  in 
this  configuration  is  perpendicular  to  the  exhaust  flow,  the  thermal  conduction  in  the  flow 
direction  will  be  inhibited.  Note  that  the  inflow  speed  Vjn,  which  is  a  factor  that  governs  the 
CS  geometry,  is  calculated  from  rVinBcor  =constant  under  the  steady  state  assumption  and 
with  vin  =  O.OIua  at  the  DR. 

Figure  4  plots  the  resulting  ne(r),  Te(r)  and  vout(r)  for  the  3  cases  of  the  coronal  ne/Te 


profiles,  and  for  DR  at  1.1  and  1.5  R0  for  each  case.  The  profiles  from  Vrsnak  et  al.  (2009) 
(‘lMK_Parker’)  and  the  electron  densities  derived  from  UVCS  and  LASCO  data  are  also 
plotted  for  comparison.  The  profiles  with  DR  heights  in  between  1.1  and  1.5  R0  lie  between 
the  two  curves  for  each  case  respectively.  Note  that  models  N1T1  and  N1T2  have  almost  the 
same  profiles  in  the  CS  because  when  beta  is  small  as  in  these  cases,  the  jump  condition  for 
ne  is  about  the  same  (Eq.(l)).  The  factor  of  2  difference  in  /3  is  compensated  by  the  factor 
of  2  in  Tj  (Eq.(2)  implies  Ti/T\  ~  /3_1  for  small  f3),  leading  to  almost  the  same  temperature 
T2  in  the  CS.  The  electron  density  in  the  CS  is  higher  than  that  in  the  ambient  corona  at 
all  heights  (comparing  Figs. 3  and  4,  see  also  Eq.(l)),  therefore  the  CS  structure  can  stand 
out  against  the  ambient  corona  as  observed. 


In  order  to  calculate  the  emission  from  the  CS,  we  need  to  first  calculate  the  ion  charge 
states  (i.e.  ionic  fractions)  of  the  ions  that  contribute  to  the  emission.  With  large  outflows 
in  the  CS,  the  ion  charge  states  do  not  always  maintain  ionization  equilibrium  when  the 
electron  density  is  low  enough.  The  evolution  of  the  ion  charge  states  with  the  flow  depends 
on  the  electron  density  and  temperature  profiles  as  well  as  the  ion  outflow  profiles.  For  a 
given  element,  the  evolution  of  ion  charge  states  with  the  flow  is  expressed  as  (under  steady 
state  assumption): 


dyq 

dt 


lVq—  1  C'i^qVq)  “I-  (-Rrr,g+ 1  Rdr,q-\- 1)  Vq-\- 1  (Rrr,q  Rdr,q)  Vq 


(10) 


where  r  is  the  heliocentric  height  of  the  fluid  element  along  the  CS,  yq  is  the  ionic  fraction 
of  charge  state  q.  Ci:q,  Rrr,q,  Rdr,q  (which  mainly  depend  on  Te )  are  the  rates  for  electron 
impact  ionization  (including  auto-ionization),  radiative  recombination  and  dielectronic  re¬ 
combination  respectively,  out  of  the  charge  state  q.  The  ionization  and  recombination  rates 
are  calculated  using  the  most  recent  compilation  by  Bryans  et  al.  (2006,2009  and  references 
within).  Here  we  assume  that  all  ions  have  the  same  outflow  speed.  For  each  element, 
Eq.(10)  for  all  charge  states  (i.e.  Z+l  equations  for  an  element  with  atomic  number  Z)  are 
solved  simultaneously.  We  calculate  the  ionic  fractions  for  13  most  abundant  elements:  H, 
He,  C,  N,  O,  Ne,  Mg,  Si,  S,  Ar,  Ca,  Fe,  Ni.  Note  that  the  ions  that  just  cross  the  SMS  should 
still  carry  the  charge  state  distribution  at  the  ambient  coronal  temperature  (even  though  the 
electron  temperature  jumps  to  T2  after  the  crossing),  and  those  calculated  by  Eq.(10)  are  for 
ions  coming  from  the  cell  below  (cp.  Fig.lc)  following  the  flow.  Therefore,  the  charge  state 
distribution  in  a  given  cell  would  be  the  average  of  the  two  regions  weighted  by  the  mass 
inside  Vsms  and  Vthru ,  respectively,  for  that  given  cell. 


Figure  5  plots  the  resulting  yq{r)  for  case  N1T2  with  DR  at  1.2  R0  for  a  few  ions, 
along  with  the  values  in  the  case  of  ionization  equilibrium.  We  can  see  that  these  ion  charge 
states  are  far  from  ionization  equilibrium  above  a  certain  height  (different  for  different  ion 
species),  and  the  ionic  fractions  ‘freeze- in’  as  they  do  in  the  solar  wind.  Therefore,  to  predict 
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emissions  inside  the  CS,  it  is  important  to  take  this  non-equilibrium  condition  into  account, 
and  it  is  not  necessarily  valid  to  calculate  the  emission  based  on  ionization  equilibrium  (i.e. 
only  based  on  the  electron  temperature).  Similarly,  the  electron  temperature  derived  from 
line  ratios  assuming  ionization  equilibrium  may  not  be  the  actual  electron  temperature. 

Figure  6  compares  yq(r)  for  the  three  cases  at  DR  height  of  1.2  R0  for  a  few  ions.  We 
can  see  that  there  are  significant  differences  between  models  of  different  coronal  profiles  that 
would  result  in  different  emission  intensities  at  various  heights,  as  well  as  different  frozen-in 
charge  states.  Note  that  even  though  models  N1T1  and  N1T2  have  almost  the  same  ne ,  Te, 
and  vout  (Fig.4),  the  different  ambient  coronal  temperatures  result  in  different  initial  charge 
states  at  the  location  of  the  SMS  crossing,  thus  different  evolution  of  the  ionic  fractions  as 
the  ions  flow  along  the  CS.  Therefore  the  observations  can,  in  principle,  be  used  to  constrain 
these  input  model  parameters.  Note  that  many  ions  seem  to  start  freezing-in  at  higher 
heights  than  what  are  usually  expected  for  the  normal  solar  wind  (e.g.  Biirgi  &  Geiss  1986, 
Ko  et  all.  1997),  probably  due  to  much  higher  density  inside  the  CS  than  in  the  solar  wind. 


2.2.  Model  2:  The  Streamline  Model 


Another  model  we  explored  assumes  that  the  plasma  that  just  comes  inside  the  CS 
through  SMS  forms  its  own  ‘streamlines’  and  does  not  mix  with  the  material  that  comes  in 
at  other  SMS  crossings.  The  properties  inside  the  SMS  at  a  given  location  would  then  be  an 
average  quantity  from  these  streamlines  weighted  by  the  volume  each  of  these  streamlines 
occupies  at  that  location  (see  below).  While  the  long  collisional  mean  free  path  suggests  the 
idea  of  mixing  as  in  Model  1,  it  is  not  clear  if  the  actual  magnetic  field  configuration  inside 
the  SMS  is  the  same  as  that  in  the  Petchek’s  model.  Also,  MHD  turbulence  may  inhibit  the 
mixing  of  the  plasma  and  thermal  equilibration  (turbulent  speeds  of  ~60  krn/s  have  been 
found  in  the  post-CME  CSs  by  Bemporad  (2008)).  Therefore  we  also  want  to  examine  how 
such  a  ‘Streamline  model’,  which  can  be  regarded  as  an  extreme  opposite  case  of  Model  1, 
compares  with  the  observations. 


We  use  the  same  three  coronal  models  (i.e.  N1T1,  N1T2,  N2T2),  and  the  same  com¬ 
pression/heating/acceleration  at  the  SMS  crossing  (Eqs.(l),  (2),  (4)).  As  the  plasma  flows 
in  each  streamline  after  the  SMS  crossing  (cp.  Fig. lb),  the  electron  density  decreases  due  to 
spherical  expansion: 


ne{r)  =  nejSMS 


r2 

'  SMS 


(11) 


where  n£tsMS  (same  as  n2  in  Eq.(l))  is  the  electron  density  at  the  SMS  crossing  at  height 
r sms-  The  evolution  of  the  electron  temperature  is  governed  by  adiabatic  cooling  and  ra- 


diative  cooling: 


(12) 


dTejr)  _  _4 Tg  _  2nenHF 
dr  3  r  3(1.92 neku(r))' 

where  A  is  the  radiative  cooling  rate.  Because  the  ions  are  most  likely  not  in  ionization 
equilibrium  along  the  flow  (cp.Eq.(lO)  and  Fig.5),  we  use  the  non-equilibrium  ionic  fractions 
(Eq.(10))  to  calculate  the  radiative  cooling  rate.  The  radiative  cooling  includes  bound-bound, 
bound-free  and  free-free  processes  using  routines  provided  by  CHIANTI  atomic  database  ver¬ 
sion  6.0  (Dere  et  al.  1997,  2009)  but  with  ionic  fractions  calculated  here  replacing  those  in 
ionization  equilibrium.  We  find  that  the  radiative  cooling  is  negligible  compared  to  adia¬ 
batic  cooling  in  almost  all  cases  except  at  the  few  lowest  heights  where  the  density  is  high. 
Therefore  it  can  be  neglected  in  general.  The  flow  speed  in  each  streamline  is  equal  to  the 
coronal  Alfven  speed  at  the  SMS  crossing  (Eq.(4)).  The  evolution  of  the  ionic  fractions  in 
each  streamline  thus  can  be  calculated  according  to  Eq.(10).  Note  that  we  neglect  heat 
conduction  in  Eq.(12)  assuming  that  MHD  turbulence  inside  the  CS  suppresses  the  thermal 
conduction.  Long  current  sheets  are  expected  to  undergo  tearing  instability,  which  can  cre¬ 
ate  a  large  number  of  plasmoids  at  different  scales  (e.g.  Shibata  &  Tanuma  2001,  Barta  et 
al.  2010,  and  references  therein).  Due  to  the  poloidal  field  of  the  plasmoids,  they  are  ther¬ 
mally  ‘insulated’  in  such  medium  that  is  filled  with  large  number  of  stochastically  moving 
plasmoids,  thus  reducing  the  thermal  conduction.  There  will  be  thermal  conduction  at  some 
level  in  spite  of  the  effects  of  turbulence,  but  we  have  chosen  the  extreme  case  to  investigate 
the  range  of  possible  results. 

Once  the  electron  density,  temperature  and  ion  outflow  profiles  in  each  streamline  are 
known,  the  ionic  fraction  profiles  in  the  streamline  can  be  calculated  according  to  Eq.(10). 
The  properties  inside  the  CS  at  a  certain  location,  in  the  context  of  this  model,  are  therefore 
average  quantities  taking  into  account  all  ‘streamlines’  that  form  between  the  DR  and  that 
given  location,  assuming  that  the  streamline  structure  inside  the  CS  can  not  be  resolved. 
Under  spherical  geometry,  the  volume  at  height  r  occupied  by  a  streamline  that  entered  the 
CS  (i.e.  crossing  the  SMS)  at  height  tsms  is  proportional  to  r/rsMS ■  The  average  quantity 
of  a  parameter  X  observed  at  height  r  is  thus 


X(r) 


^-‘rSSMS=rDRXSMs{r )F(r,  T SMS ) 


KSsTs=rDRF(r,  r sms) 


(13) 


where  tsms  and  tdr  are  the  height  of  the  SMS  crossing  and  diffusion  region,  respectively. 
Xsms  is  the  physical  quantity  at  height  r  in  the  streamline  that  enters  the  SMS  at  tsms, 
and  F(r,rsMS )  is  the  ‘contribution  factor’  from  the  streamline  forming  at  tsms ■  F(r,  tsms) 
is  different  for  different  physical  parameters.  For  the  electron  density,  F(r,rsMs)  =  f/f’SMS- 
For  the  ionic  fraction  yq ,  F(r,rsMs)  =  rne(r)/rSMS  (ion  density  is  proportional  to  yq  *  ne). 
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Figure  7  shows  an  example  of  how  different  streamlines  contribute  to  an  ionic  fraction 
observed  at  height  r  =  1.7  R0  for  the  case  of  tdr  =  1.2  R0.  All  3  coronal  models  are  plotted 
for  comparison.  Also  plotted  are  those  for  the  electron  density.  We  choose  those  ions  that 
emit  high-temperature  lines  observable  by  UVCS  in  the  CS  at  this  height  (i.e.  Fe+1'([Fe 
XVIII]  A974),  Si+11(Si  XII  A499),  Ca+13([Ca  XIV]  A943)  (e.g.  Ciaravella  et  al.  2002,  Ko 
et  al.  2003,  Bemporad  et  al.  2006).  One  can  see  that  the  streamlines  (i.e.  r sms's)  that 
make  a  major  contribution  to  the  average  ionic  fraction  are  different  for  different  ions  and 
different  coronal  models,  but  most  of  the  contribution  is  from  SMS  crossings  at  low  heights. 
Similar  examination  for  an  observed  height  of  1  AU  indicates  that  the  major  contribution 
for  most  abundant  ions  at  1  AU  is  from  streamlines  that  enter  the  SMS  below  2  R0.  For  the 
electron  density,  it  is  the  same  for  all  3  models  because  ne(r )  evolves  in  the  same  way  with 
r  (Eq.(ll)). 

Figure  8  plots  the  average  ionic  fractions  (Eq.(13))  for  selected  ions  as  a  function  of 
height  for  the  3  coronal  models.  Comparing  with  the  Fully-mixed  model,  we  find  that  the 
results  of  the  two  models  are  similar  within  an  order  of  magnitude.  Figure  9  plots  the 
average  electron  density  compared  with  the  data  (cp.Fig.4).  Fig.9  shows  that,  similar  to  the 
Fully-mixed  model,  models  N1T1  and  N1T2  with  reasonable  DR  heights  predict  electron 
densities  that  agree  with  the  observations,  while  model  N2T2  seems  to  underestimate  the 
electron  density. 


3.  Predicted  Observables  from  Model  Calculations 

The  main  purpose  of  this  work  is  to  test  post-CME  CS  models  and  see  if  model  pre¬ 
dictions  agree  with  observations.  Figs. 4  and  9  show  that  the  predicted  electron  density 
profiles  within  the  CS  overall  agree  with  the  data  for  both  model  approaches  (see  also  Fig. 7 
of  Vrsnak  et  al.  2009).  The  two  models  predict  similar  electron  densities  at  low  heights,  but 
the  electron  density  profile  falls  faster  with  height  for  the  Fully-mixed  model  at  r  ^  4  R0. 
Therefore  CS  observations  over  an  extended  range  of  heights  may  be  able  to  distinguish  such 
different  model  approaches  and  indicate  which  one  represents  better  the  actual  conditions 
inside  the  CS.  It  is  also  possible,  in  the  context  of  either  model,  to  constrain  the  coronal 
parameters  outside  of  the  CS. 
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3.1.  UV  Line  Emission 

Another  predicted  quantity  is  the  line  emission  (see  Sec.l).  In  coronal  conditions,  the 
electron  density  is  low  enough  that  most  emission  lines  are  produced  by  electron  collisional 
excitation  followed  by  spontaneous  emission.  For  a  plasma  with  electron  temperature  Te, 
the  line  emission  is  thus: 

hue  =  7-—  [  G(Te)dEM(Te)  photon  cm-2s-1sr-1  (14) 

4 vr  nH  J 

where  nei/riH  is  the  elemental  abundance  relative  to  hydrogen  (i.e.  absolute  abundance, 
Grevesse  et  al.  2007).  G(Te )  is  the  contribution  function  which  is  defined  as 

77  • 

G(Te)  =  —Blineqline(Te)  (15) 

riel 

where  num/nei  is  the  ionic  fraction  calculated  from  Eq.(10).  Bune  is  the  branching  ratio  for 
the  line  transition,  and  qune(Te )  is  the  electron  excitation  rate  which  is  a  function  of  only 
the  electron  temperature  in  the  low-density  limit.  dEM(Te )  =  d(neriHL )  is  the  emission 
measure  (in  cm-5)  at  a  given  electron  temperature  with  line-of-sight  (LOS)  depth  L,  and 
Uh  —  ne/1.2  for  fully  ionized  plasma.  In  our  models,  the  ionic  fractions  are  not  in  ionization 
equilibrium  (e.g.  see  Fig.5).  Therefore  we  calculate  G  from  the  CHIANTI  atomic  database 
version  6.0  (Dere  et  al.  1997,  2009)  but  with  the  ionic  fractions  calculated  here  replacing 
those  in  ionization  equilibrium  (e.g.  Figs. 6  and  8). 

For  the  Fully-mixed  model,  the  line  emission  can  be  calculated  directly  from  Eq.(14) 
with  the  electron  density,  temperature  and  ionic  fractions  calculated  in  the  model.  For  the 
Streamline  model,  one  needs  to  calculate  the  average  quantity  as  in  Eq.(13)  with  XsMs(r)  = 
Iiine,SMs(r)  and  F(r,rSMS )  =  r/rSMS ■  There  are  two  assumptions  needed.  One  is  the 
elemental  abundance,  and  the  other  is  the  LOS  depth.  The  First  Ionization  Potential  (FIP) 
effect  is  commonly  observed  in  coronal  plasmas.  The  FIP  bias,  defined  as  the  abundance 
ratio  of  a  low-FIP  element  (FIP  <  10  eV)  to  a  high-FIP  element  (FIP  >  10  eV)  relative  to 
its  photospheric  ratio,  is  generally  in  the  range  of  3-4  in  active  regions  and  streamers  (e.g. 
see  review  by  Raymond  et  al.  2001  and  references  within).  Since  fast  CMEs  with  such  WL 
ray/CS  structure  almost  always  occur  from  active  regions,  we  assume  that  the  abundance 
for  all  low-FIP  elements  is  3x  its  photospheric  value,  and  the  abundance  for  all  high-FIP 
elements  is  equal  to  its  photospheric  value  (Grevesse  et  al.  2007).  For  the  LOS  depth,  we 
adopt  the  value  of  0.05  R0  at  1.1  R0  so  the  LOS  depth  at  a  given  height  r  is  0.05  x  r/1.1 
R0  under  spherical  expansion.  For  individual  CME  events,  this  LOS  depth  can  be  estimated 
from  the  length  and  inclination  of  the  neutral  line  at  the  eruption  site  (Ciaravella  &  Raymond 
2008,  Vrsnak  et  al.  2009,  Lin  et  al.  2009).  In  any  case,  different  assumptions  for  the  FIP 
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bias  and  LOS  depth  can  be  easily  taken  into  account  since  the  calculated  line  fluxes  are 
linearly  proportional  to  these  two  quantities. 

Figure  10  plots  the  predicted  fluxes  for  Si  XII  A499,  [Ca  XIV]  A943  and  [Fe  XVIII]  A974 
versus  heliocentric  heights  for  the  Fully-mixed  model  with  3  coronal  models  and  DR  height 
at  1.2  R.0.  As  previously  mentioned,  these  lines  are  chosen  because  they  were  commonly 
observed  by  LIVCS  in  post-CME  CS  events.  UVCS  data  for  4  events  are  also  plotted  for 
comparison.  Fig. 11  plots  the  same  but  for  the  Streamline  model.  The  two  models  predict 
similar  line  fluxes  at  low  heights  (  A-  2.0  R0).  The  difference  is  at  higher  heights  where  the 
line  fluxes  in  the  Fully-mixed  model  decrease  more  rapidly  with  height.  At  5  R0,  these  lines 
are  fainter  by  an  order  of  magnitude  in  the  Fully-mixed  model  than  the  Streamline  model. 
At  heights  of  the  available  UVCS  observations  (<  1.7  R0),  the  model  predictions  agree  with 
the  observations  in  an  overall  sense.  The  Mar. 23,  1998  event  seems  to  be  more  consistent 
with  the  N1T1/N1T2  models,  and  the  Jan. 8,  2002  event  seems  to  be  more  consistent  with  the 
N2T2  model.  This  would  indicate  different  ambient  coronal  conditions  for  different  events. 
For  the  Nov.4,  2003  event,  however,  Si  XII  and  Fe  XVIII  lines  are  consistent  with  different 
coronal  models.  One  possible  explanation  is  that  none  of  the  coronal  models  adopted  here 
represents  the  actual  coronal  condition  for  this  event.  A  higher  temperature  in  CS  (e.g. 
from  smaller  coronal  /3,  Eq.(2))  would  produce  less  Si  XII  emission  relative  to  the  Fe  XVIII 
emission.  For  Ca  XIV,  all  models  underestimate  the  line  emission  up  to  a  factor  of  5.  This 
line  is  blended  with  [Si  VIII]  A944  but  the  absence  of  low-temperature  emission  along  the 
slit  where  the  Fe  XVIII  emission  exists  rules  out  line  blending  as  a  possible  cause.  The 
discrepancy  may  lie  in  the  atomic  data  and  the  assumed  coronal  abundance. 

Figure  12  shows  an  example  of  how  different  streamlines  contribute  to  the  line  fluxes 
observed  at  height  r  =  1.7  R0  for  the  case  of  tdr  =  1.2  R0.  All  3  coronal  models  are 
plotted  for  comparison.  Depending  on  the  individual  line,  the  dominant  emission  comes 
from  ‘streamlines’  of  a  limited  range  of  the  SMS  crossings.  There  is  little  difference  among 
the  3  coronal  models.  This  is  in  contrast  with  the  ionic  fractions  (Fig. 7),  probably  because 
of  a  combination  of  the  effects  of  the  electron  density  and  temperature  on  the  line  emission. 

One  interesting  aspect  to  examine  is  the  change  of  line  emission  with  time.  As  the 
reconnection  continues,  the  location  of  the  DR  will  move  up  toward  higher  height  with 
speed  approximately  equal  to  the  inflow  speed  of  the  reconnection  (Lin  &  Forbes  2000). 
Since  the  formation  of  the  SMS  is  much  faster  (in  Alfven  speed)  than  the  movement  of  the 
DR,  the  SMS  developed  from  a  given  DR  height  can  be  regarded  as  a  steady  state  structure. 
The  properties  within  the  CS  will  then  change  with  the  movement  of  the  DR.  The  time 
evolution  of  these  properties  may  provide  information  about  the  DR  movement,  therefore 
the  inflow  speed/reconnection  rate.  Figures  13  and  14  plot  the  3  line  fluxes  and  the  electron 
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density  versus  DR  height  as  seen  at  3  locations  (1.5,  1.6,  1.7  R0),  along  with  the  UVCS 
data,  for  the  Fully-mixed  model  and  Streamline  model,  respectively.  The  time  sequence  of 
line  fluxes  and  electron  density  in  the  Nov. 4,  2003  event  (purple  symbols,  Ciaravclla  et  al. 
2008,  observed  at  1.66  R0)  indicates  that  during  2.5  hours,  the  Fe  XV111  fluxes  decreased  by 
a  factor  of  0.63,  Si  Xll  fluxes  increased  by  a  factor  of  1.5,  and  the  electron  density  decreased 
by  a  factor  of  0.7.  If  we  use  model  N1T2  as  a  proxy  (note  that  the  direction  of  change  in 
these  3  quantities  is  the  same  as  the  data,  even  though  the  magnitude  does  not  all  fit),  this 
roughly  corresponds  to  an  upward  motion  of  around  0.15  R0  which  in  turn  corresponds  to 
an  inflow  speed  of  ~10  krn/s  and  M 4  ~0.01  (based  on  Alfven  speed  of  ~  O(103)  krn/s). 
This  is  consistent  with  both  theoretical  (although  on  the  low  side,  e.g.  Lin  &  Forbes  2000) 
and  observational  (e.g.  Yokoyama  et  al.  2002,  Ko  et  al.  2003,  Lin  et  al.  2005)  findings. 


3.2.  Simulated  Broad-Band  X-Ray  Emission 

Besides  particular  lines  of  interest  shown  above,  our  model  of  the  electron  density, 
temperature  and  ionic  fractions  for  the  13  most  abundant  elements  enables  us  to  calculate 
emission  spectra  across  a  continuous  wavelength  range.  One  such  application  is  to  use  the 
emission  code  of  Raymond  &  Smith  (1977;  updated  in  Cox  &  Raymond  1985)  which  calcu¬ 
lates  the  line  and  continuum  emission  spectrum.  This  predicted  spectrum  can  then  be  com¬ 
pared  with  observations  from  broad-band  instruments  such  as  Hinode/XRT  and  SDO/AIA, 
after  folding  with  the  instrument’s  wavelength  response  function,  for  specific  CS  events  these 
instruments  might  observe. 

To  demonstrate  this,  we  use  the  non- equilibrium  ionization  states  computed  with  both 
models  along  with  the  elemental  abundances  described  above  to  compute  the  emission  over 
the  range  2.5-4998  eV  (2.48-4960  A).  The  Raymond  &  Smith  model  includes  fewer  emission 
lines  than  the  APEC  code  (Smith  et  al.  2001),  but  for  low  spectral  resolution  data  the 
results  are  very  similar.  We  then  fold  the  emission  spectrum  with  the  wavelength  response 
function  of  the  XRT  instrument  on  Hinode  to  obtain  the  predicted  count  rate  of  the  CS 
if  it  were  to  be  seen  by  XRT.  The  wavelength  response  function  is  calculated  using  the 
standard  SolarSoft  routines  with  contamination  thickness  for  Apr. 9,  2008,  14  LIT.  We  chose 
this  particular  instrument  and  date  because  a  very  interesting  CME/CS  event  was  observed 
by  SOHO  and  Hinode  on  this  day  that  has  been  analyzed  in  great  detail  (Savage  et  al.  2010, 
Landi  et  al.  2010). 

Figures  15  plots  the  predicted  count  rate  (DN/sec/pixcl)  versus  height  for  the  Fully- 
mixed  model  as  seen  by  Hinode/XRT  thin  Al-Poly  filter  for  the  3  coronal  models  and  3  DR 
heights.  Figure  16  plots  the  same  for  the  Streamline  model.  The  LOS  depth  is  assumed  to  be 


15 


the  thickness  of  the  CS,  5  x  103  km,  as  determined  by  Savage  et  al.  (2010)  for  the  Apr. 9,  2008 
event.  We  can  see  that  the  count  rate  can  differ  by  an  order  of  magnitude  among  coronal 
models  and  is  mainly  dictated  by  the  ambient  coronal  density,  not  the  electron  temperature. 
Figure  17  shows  the  CS  observed  by  Hinode/XRT  on  Apr. 9,  2008.  Plotted  in  the  upper 
panels  are  the  measured  count  rates  along  the  solar  Y-position  at  4  solar  X-positions  that 
cross  the  CS  (marked  in  the  lower  panels).  The  location  of  the  DR  is  at  ~1.25  R0  based 
on  the  downflows/upflows  measured  by  Savage  et  al.  (2010).  This  is  around  the  second 
outermost  X-position  marked  in  Fig.  17  (inner  solid  line).  The  ‘bumps’  between  Y  pixel  of 
100-300  are  where  the  CS  is  located.  Plotted  in  red  on  the  two  lowest  curves  (solid  lines)  are 
the  Y-pix  ranges  chosen  to  estimate  the  mean  and  standard  deviation  of  the  count  rate  within 
the  CS.  These  numbers  are  shown  as  blue  data  points  on  Figs. 15  and  16.  The  heliocentric 
heights  of  these  data  points  were  estimated  based  on  that  the  flare/CME  source  was  ~  23° 
behind  the  limb  at  14  UT  on  Apr. 9,  2008  (Savage  et  al.  2010).  Comparing  the  predicted 
and  observed  count  rates  at  these  X-positions  indicates  that  our  model  calculations  are  in 
the  agreeable  range  with  the  data  for  the  N1T1  and  N1T2  models,  but  not  for  the  N2T2 
model.  Note  that  our  predicted  count  rates,  bracketing  the  data  between  DR  height  of  1.2 
and  1.3  R0,  are  also  consistent  with  the  position  of  the  X  point  (~  1.25  R0)  as  found  out 
by  Savage  et  al.  (2010).  The  agreement  implies  that  the  ambient  coronal  density  profile  is 
close  to  the  ‘Nl’  profile,  at  least  at  around  1.3  R0.  It  should  be  possible  to  compare  the  ‘Nl’ 
profile  with  the  electron  density  profile  derived  from  the  white  light  polarization  brightness 
(pB)  measurement  at  the  time  before  and  during  the  event  and  see  if  they  agree,  but  we  do 
not  make  such  attempt  in  this  paper. 


3.3.  Solar  Wind  Charge  States  at  1  AU 

Many  of  the  magnetic  clouds  (MCs)  seen  as  interplanetary  coronal  mass  ejections 
(ICMEs)  bear  high  ion  charge  states  usually  represented  in  high  0+7/0+6  ratios  (>  0.2) 
and  high- ionization  charge  states  of  Fe,  e.g.  Fe+16  (Lepri  &  Zurbuchen  2004).  Since  the 
ion  charge  states  usually  freeze-in  near  the  Sun,  this  implies  electron  heating  in  the  CME 
material  in  the  early  stage  of  the  eruption  process  (e.g.  Akmal  et  al.  2001,  Rakowski  et 
al.  2007,  Lee  et  al.  2009).  Our  model  of  the  post-CME  CS  also  predicts  higher  charge 
states  due  to  higher  electron  temperature  than  the  ‘normal’  corona  from  the  electron  heat¬ 
ing  crossing  the  SMS,  as  shown  in  our  results.  Plasma  ejected  upwards  in  the  current  sheet 
enters  the  flux  rope  formed  by  reconnection  below  the  CS  core.  If  there  is  a  pre-existing 
flux  rope,  reconnection  forms  a  flux  rope  around  it  having  a  comparable  amount  of  magnetic 
flux,  and  therefore  volume  (Lin  et  al.  2004,  Mostl  et  al.  2008),  so  that  much  of  the  plasma 
in  a  magnetic  cloud  has  passed  through  the  current  sheet.  However,  there  is  no  report  so 
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far  of  definite  signatures  of  post-CME  CS  measured  in-situ  at  1  AU.  Note  that  this  ‘CS’  is 
different  from  the  ‘in-situ  reconnection  exhaust’  reported  by  Gosling  et  al.  (2005),  but  the 
model  work  presented  here  can  in  principle  be  applied  to  this  case  with  very  different  in-situ 
ambient  conditions  outside  the  SMS. 

Figures  18  and  19  show  the  predicted  frozen-in  0+7/0+6  ratios  and  the  average  Fe 
charge  at  1  AU  for  the  Fully-mixed  model  and  Streamline  model,  respectively.  We  can  see 
that  there  are  notable  differences  between  the  two  model  approaches,  as  well  as  the  3  coronal 
models.  In  general,  the  Streamline  model  predicts  higher  charge  states  at  1  AU  (or  >  20 
R0  according  to  Fig.6  and  8)  than  the  Fully-mixed  model.  The  ionization  charge  state  is 
higher  for  model  N2T2  than  the  other  two  coronal  models  due  to  higher  CS  temperature 
(e.g.  Fig.4),  as  opposed  to  that  the  emission  is  lower  for  model  N2T2  due  to  lower  density. 
We  note  that  the  range  of  these  predicted  quantities  of  0+7/0+6  and  average  Fe  charge 
is  common  in  the  ICME  solar  wind  data  (0+'/0+6  >  0.2,  average  Fe  charge  >  12.  See, 
e.g.  Lepri  &  Zurbuchen  2004,  Wimmer-Schweingruber  et  al.  2006,  Zurbuchen  &  Richardson 
2006).  If  some  of  the  ICMEs  with  high  0+7/0+6  and  Fe  charge  states  are  associated  with 
such  post-CME  CS  structure,  It  would  imply  that  the  diffusion  region  remains  low  in  the 
corona.  It  remains  to  be  seen  if  the  models  presented  here  can  be  tested  by  observations  of 
charge  states  within  the  to-be-identified  post-CME  CS  structure  in  the  solar  wind. 


4.  Discussion 

In  order  to  understand  post-CME  WL  rays  and  current  sheets,  it  is  important  to  have 
quantitative  comparisons  between  theory  and  observations.  The  models  presented  here  are 
a  first  step  in  that  direction.  Several  crucial  observations  exist,  including  the  Fe  XXIV 
spines  seen  above  some  post-flare  loops  in  TRACE  images  (Innes  et  al.  2003a),  the  high 
velocity  Fe  XIX  and  Fe  XXI  lines  seen  by  SUMER  in  several  events  (Innes  et  al.  2003b; 
Wang  et  al.  2007),  the  UVCS  observations  of  highly  ionized  gas  (Ciaravella  et  al.  2002,  Ko 
et  al.  2003,  Bemporad  et  al.  2006,  Ciaravella  &  Raymond  2008),  evidence  of  reconnection 
downflows/upflows  observed  by  Hinode/XRT  (Savage  et  al.  2010),  and  in  situ  measurements 
of  high  ionization  states  in  ICMEs  (e.g.  Lepri  &  Zurbuchen  2004). 

We  have  taken  a  simple  model  of  the  exhaust  region  of  Petschek  reconnection  and 
computed  the  time-dependent  ionization  state  in  order  to  predict  UV  emission  line  intensities 
and  the  count  rates  in  X-ray  bands.  This  model  is  consistent  with  observations,  thus  it  is 
a  viable  explanation  for  the  observed  ‘post-CME  current  sheet’.  This  does  not  exclude 
other  possible  models  but  any  model  should  be  able  to  predict  properties  that  can  produce 
observables  that  are  consistent  with  the  observations.  The  Vrsnak  et  al.  (2009)  model,  in 
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our  opinion,  is  the  most  viable  model  so  far  that  adequately  treats  plasma  properties  of 
the  CS  above  the  DR  for  which  we  can  calculate  physical  parameters  and  compare  with  the 
observations.  Given  the  model  assumptions,  the  good  agreement  with  the  observations  in 
both  density  and  UV  line /X-ray  emission  is  more  than  just  fortuitous.  In  the  context  of  this 
model,  the  width  of  the  WL  ray  increases  with  height  which  is  also  consistent  with  the  data 
(see  Vrsnak  et  al.  2009)  .  This  indicates  that  this  model  can  be  the  right  interpretation 
of  the  observed  WL  ray/CS.  In  principle,  our  calculations  can  be  applied  to  any  post-CME 
CS  models  that  provide  information  about  the  electron  density/temperature  and  ion  flow 
speed,  such  as  in  the  MHD  models  (e.g.  Riley  et  al.  2008  by  using  proton  parameters).  Note 
that  the  main  purpose  of  this  paper  is  to  show  how  this  particular  post-CME  CS  model 
concept  agrees  with  the  observations  but  we  do  not  attempt  to  make  fine  adjustments  of  the 
model  parameters  to  fit  data  of  any  one  particular  CS  event.  It  is  also  possible  to  apply  such 
calculations  to  model  the  downflow  region  below  the  DR  (e.g.  Yokoyama  &  Shibata  1998, 
2001)  but  we  do  not  make  such  attempt  in  this  paper. 

Both  CS  models  (‘Fully-mixed’  and  ‘Streamline’)  are  found  to  give  similar  UV  and  X-ray 
emission  at  heights  below  ~2  R0  but  it  drops  more  quickly  with  height  in  the  Fully-mixed 
model.  Thus  it  seems  that  CS  observations  higher  than  2  R0  would  be  able  to  differentiate 
the  two  models.  However,  current  instruments  do  not  have  the  sensitivity  yet  to  detect  such 
high  temperature  emission  much  beyond  what  is  presented  here.  Figs. 4  and  9  indicate  that 
electron  densities  derived  from  WL  observations  beyond  4  R0  can  be  used  to  differentiate 
the  two  models.  Such  data  are  available  although  we  do  not  attempt  to  acquire  in  this  work. 
The  ionic  charge  state  (Figs. 18  and  19)  is  another  useful  parameter  to  differentiate  these 
models.  However,  an  unambiguous,  positive  identification  of  such  post-CME  CS  in  the  solar 
wind  is  necessary  before  such  comparison  can  be  made. 

The  model  results  depend  strongly  on  the  model  of  the  ambient  corona.  The  differences 
among  the  three  coronal  models  show  that,  within  the  context  of  the  models  here,  the 
coronal  electron  density  profile  is  the  most  important  factor  in  determining  the  emission 
properties  inside  the  CS.  In  this  work,  we  did  not  investigate  the  effect  of  the  ambient 
coronal  magnetic  held.  But  we  expect  it  will  also  play  a  crucial  role  (as  dictating  the  plasma 
P)  in  determining  the  CS  properties.  The  more  accurately  the  ambient  conditions  at  the 
time  of  the  CS  event  are  known,  the  better  such  observations  can  pinpoint  the  reconnection 
processes  and  properties  in  the  post-CME  CS.  And  comparison  of  such  models  presented 
here  with  observations  can  offer  a  possible  diagnostic  tool  for  inferring  the  external  coronal 
conditions  of  the  post-CME  CS. 

As  a  final  remark,  the  calculation  of  plasma  emission  and  time-dependent  ionization 
rely  on  atomic  rates  for  collisional  excitation,  ionization  and  recombination.  Aside  from  the 
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accuracy  issue  of  these  atomic  rates,  the  electron  velocity  distribution  function  (VDF)  is  an 
important  factor  in  these  rate  calculations.  In  this  work,  we  assume  that  the  electron  VDF  is 
a  Maxwellian.  The  deviation  of  the  VDF  from  a  Maxwellian  can  affect  these  rates  in  different 
degrees  for  different  ions  (e.g.  Ko  et  al.  1996,  Dzifcakova  &  Mason  2008).  One  future  effort 
of  our  calculations  will  be  to  investigate  possible  effects  of  a  non-Mawellian  electron  VDF.  A 
kinetic  approach  in  magnetic  reconnection/CS  models  that  can  provide  information  about 
the  electron  VDF  would  be  helpful  in  nailing  down  our  knowledge  of  the  post-CME  CS  with 
observations. 
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Fig.  1. —  (a)  Standard  flare-CME  model  depicting  a  current  sheet  between  the  post-flare 
loops  and  the  ejecting  flux  rope  (e.g.  Lin  &  Forbes  2000).  (b)  Sketch  of  the  model  by  Vrsnak 
et  ah  (2009)  that  interprets  the  observed  WL  ray  (CS)  as  the  feature  bounded  by  a  pair 
of  slow-mode  shocks  (SMS)  above  the  diffusion  region  (DR),  (c)  Sketch  of  the  ‘Fully-mixed 
model’  (adopted  from  Vrsnak  et  al.  2009). 
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Fig.  2. —  Upper  panels:  the  CME/CS  event  on  Jan. 08,  2002  (see  Ko  et  al.  2003).  Lower 
panels:  the  CME/CS  event  on  Nov.04,  2003  (see  Ciaravclla  &  Raymond  2008).  Images 
are  from  SOHO/LASCO  (coronal  images  in  red-orange),  SOHO/EIT  A195  (solar  images  in 
green),  SOHO/UVCS  (slit  images  marked  in  the  two  right  panels),  SOHO/CDS  (coronal 
loop  image  in  grey  in  the  upper-right  panel)  and  MLSO/MK4  (coronal  image  in  blue  in  the 
upper-right  panel). 
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Fig.  3. —  Profiles  for  coronal  magnetic  field  strength  (upper  left),  electron  density  (upper 
right),  electron  temperature  (lower  left)  and  plasma  beta  (lower  right)  adopted  for  the  study. 
Also  plotted  is  the  electron  density  and  plasma  beta  profiles  for  the  ‘lMK_Parker’  model 
adopted  in  Vrsnak  et  al.  (2009).  The  electron  temperature  profile  for  the  lMK_Parker  model 
is  isothermal  at  Te(r)  =  1  x  106  K. 
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Fig.  4. —  ne(r),  Te(r)  and  vout(r)  for  the  3  cases  of  the  coronal  ne/Te  profiles,  and  for  DR 
at  1.1  (in  black)  and  1.5  R0  (in  red)  for  each  case.  The  profiles  from  Vrsnak  et  ah  (2009) 
(‘lMK_Parker’,  dash-dot-dot-dot)  and  the  electron  densities  derived  from  UVCS  and  LASCO 
data  (bine  bars  and  pluses)  are  also  plotted  for  comparison.  UVCS  data  are  from  the  Mar. 23, 
1998  event  (Ciaravella  et  al.  2002)  at  1.50  R0,  Nov. 26,  2002  event  (Bemporad  et  ah  2006) 
at  1.61  R0,  and  the  Nov. 4,  2003  event  (Ciaravella  &  Raymond  2008)  at  1.66  R0.  LASCO 
data  are  from  the  Jan. 8,  2002  event  (Ko  et  ah  2003)  at  3.0  and  4.4  R0.  Note  that  since, 
with  the  same  coronal  field,  vout  only  depends  on  the  ambient  density,  it  is  the  same  for  cases 
N1T1  and  N1T2. 
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Fig.  5. —  Ionic  fraction  profiles  for  DR  height  at  1.2  R0  for  the  Fully- mixed  model  with 
coronal  model  N1T2,  compared  with  those  in  ionization  equilibrium  at  the  local  electron 
temperature  in  the  CS  (dotted,  dash,  dash-dot  in  the  order  of  ascending  charge  state). 
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Fig.  6. —  Ionic  fraction  profiles  for  DR  height  at  1.2  R0  for  the  Fully-mixed  model.  Plotted 
are  for  3  coronal  models,  N1T1  (green-dash),  N1T2  (black-solid),  and  N2T2  (red-dash-dot). 
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Fig.  7. —  Contribution  to  the  total  ionic  fraction  observed  at  height  1.7  R0  from  streamlines 
entering  the  SMS  at  heights  from  1.2  R0  (DR  height)  to  1.7  R0.  Plotted  are  for  3  ions 
and  the  electron  density  for  3  coronal  models,  N1T1  (green-dash),  N1T2  (black-solid),  N2T2 
(red-dash-dot).  The  SMS  height  step  is  0.01  R0. 
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Fig.  8.—  Ionic  fraction  profiles  for  DR  height  at  1.2  R0  for  the  Streamline  model.  Plotted 
are  for  3  coronal  models,  N1T1  (green-dash),  N1T2  (black-solid),  N2T2  (red-dash-dot). 
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Fig.  9. —  Electron  density  profiles  for  DR  height  at  1.1  (black)  and  1.5  (red)  R0  for  the 
Streamline  model.  Plotted  are  for  3  coronal  models,  N1T1  (dash),  N1T2  (solid),  N2T2 
(dash-dot).  Also  plotted  for  comparison  are  ne  derived  from  UVCS  and  LASCO  data  (see 
caption  of  Fig.4). 
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Fig.  10. —  For  the  Fully-mixed  model:  Three  line  fluxes  versus  height  for  3  coronal  models, 
N1T1  (green-dash),  N1T2  (black-solid),  N2T2  (red-dash-dot)  for  DR  height  at  1.2  R0.  Also 
plotted  are  the  data  from  UVCS  for  the  CS  events  on  Mar. 23,  1998  (orange)  at  1.50  R0, 
Jan. 8,  2002  (blue)  at  1.53  R0,  Nov. 26,  2002  (green,  shown  is  the  range  from  Nov.27,  00:26 
UT  to  Nov. 29,  00:20  UT)  at  1.61  R0,  and  Nov. 4,  2003  (purple,  3  points  from  Nov. 4,  21:06 
UT  to  Nov.5,  00:12  UT)  at  1.66  R0.  The  coronal  abundance  is  assumed  to  be  3x  and  lx 
the  photospheric  value  for  low-FIP  and  high-FIP  elements,  respectively,  and  the  LOS  depth 
is  assumed  to  be  0.05  R0  at  1.1  R0.  Note  that  for  the  Si  XII  flux,  that  of  the  Jan. 8,  2002 
event  has  the  background  corona  taken  out.  Others  should  be  regarded  as  upper  limit. 


Si  XII  A499  Flux 


[Co  XIV]  A943  Flux 


[Fe  XVIII]  A974  Flux 


Fig.  11. —  Same  as  Fig.10  but  for  the  Streamline  model. 
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Fig.  12. —  Contribution  to  the  total  flux  observed  at  height  1.7  R0  from  streamlines  entering 
the  SMS  at  heights  from  1.2  R0  (DR  height)  to  1.7  R0.  Plotted  are  for  3  lines  for  3  coronal 
models,  N1T1  (green-dash),  N1T2  (black-solid),  N2T2  (red-dash-dot).  The  SMS  height  step 
is  0.01  R0. 


33 


Si  XII  A499  Flux  at  1.5,1. 6, 1,7  Rs 


In  10H 

E 

U 

V  1010 

U) 

c 

O 

"o 

A  io9 


10s 

1.0  1.1  1.2  1.3  1.4  1.5  1.6  1.7 

Height  of  Diffusion  Region  (Re) 


[Fe  XVIII]  A974  Flux  at  1.5,1. 6, 1.7  Rs 
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Fig.  13. —  For  the  Fully-mixed  model:  Line  fluxes  and  electron  density  observed  at  1.5,  1.6, 
1.7  Re versus  the  height  of  the  diffusion  region  for  3  coronal  models,  N1T1  (green-dash),  N1T2 
(black-solid),  N2T2  (red-dash-dot).  Also  plotted  are  the  UVCS  data  (see  Fig.  10  caption). 
For  the  3  curves  corresponding  to  the  3  observed  heights  of  each  model,  lower  height  has 
higher  flux/density  in  general  (cp.  Fig.  10). 
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Fig.  14. 


Same  as  Fig.  11  but  for  the  Streamline  model. 
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Fig.  15. —  Predicted  count  rate  vs.  height  for  Hinode/XRT  thin  Al-Poly  hlter  for  the  Fully- 
mixed  model.  Plotted  are  for  3  coronal  models,  N1T1  (green-dash),  N1T2  (black-solid), 
N2T2  (red-dash-dot),  and  for  3  DR  heights  at  1.1,  1.2,  1.3  Re.  The  LOS  depth  is  assumed 
to  be  5  x  103  km.  Also  plotted  (blue  data  points)  are  four  measurements  of  Hinode/XRT 
on  the  CS  event  on  Apr. 9,  2008  (see  Fig. 17). 
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Fig.  16. —  Predicted  count  rate  vs.  height  for  Hinode/XRT  thin  Al-Poly  filter  for  the 
Streamline  model.  See  caption  of  Fig.  15. 
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Fig.  17. —  Lower  panels:  Hinode/XRT  Al-Poly  image  on  Apr.9,  2008  at  13:46  UT  (lower 
left)  and  17:06  UT  (lower  right)  showing  the  CS  extended  above  the  post-CME  loops/cusp. 
The  vertical  lines  are  the  solar  X-positions  chosen  for  the  plots  in  the  upper  panels.  The 
spatial  binning  is  1  pixel.  Upper  panels:  observed  count  rate  in  DN/s  along  the  Y-pixel 
number  at  the  4  solar  X-positions  marked  by  the  vertical  lines  in  the  XRT  images.  The 
‘bumps’  between  Y  pixel  of  100-300  are  where  the  CS  is  located.  Note  that,  according  to 
Savage  et  al.  (2010),  the  DR  is  located  in  between  the  dashed  and  solid  lines.  Thus  our 
model  results,  which  only  apply  to  the  CS/SMS  above  the  DR,  should  be  compared  only 
with  the  data  at  the  locations  marked  by  the  two  solid  lines  (lowest  two  curves  in  the  upper 
panels).  Plotted  in  red  on  the  two  lowest  curves  are  the  Y-pixel  ranges  used  to  estimate  the 
mean  and  standard  deviation  of  the  count  rate  within  the  CS  (shown  as  blue  data  points  on 
Figs.  15  and  16). 
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Fig.  18. —  0+7/0+6  ratio  and  average  charge  of  Fe  observed  at  1  AU  for  the  Fully-mixed 
model  versus  height  of  the  DR  with  the  3  coronal  models  compared. 
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Fig.  19 


Same  as  Fig.  18  but  for  the  Streamline  model. 


